Prom Lyapunov modes to the exponents for hard disk systems 
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We demonstrate the preservation of the Lyapunov modes by the underlying tangent space dynam- 
ics of hard disks. This result is exact for the zero modes and correct to order e for the transverse and 
LP modes where e is linear in the mode number. For sufficiently large mode numbers the dynamics 
no longer preserves the mode structure. We propose a Gram-Schmidt procedure based on orthog- 
onality with respect to the centre space that determines the values of the Lyapunov exponents for 
the modes. This assumes a detailed knowledge of the modes, but from that predicts the values of 
the exponents from the modes. Thus the modes and the exponents contain the same information. 
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In a chaotic system, the difference between two nearby 
phase space trajectories, the so called Lyapunov vec- 
tor, diverges exponentially in time. If one or more of 
the rates of divergence are positive then the dynamics 
of a single initial condition is unpredictable and global 
behavior becomes important. The statistical mechan- 
ics of chaotic many particle systems is an example of a 
probabilistic treatment of the global behavior of deter- 
ministic microscopic dynamics. It is believed that the 
probabilistic axioms of statistical mechanics can be jus- 
tified by the chaotic nature of the underlying dynam- 
ics, so much effort has been devoted to finding links 
between macroscopic quantities, such as transport co- 
efficients, and chaotic properties such as the Lyapunov 
exponents P, 13, There have been some successes 
such as the conjugate pairing rule for the Lyapunov spec- 



if, III and the 



trum in some thermostated systems 
fluctuation theorem [l^, . The stepwise structure of 
the Lyapunov spectrum (the full set of Lyapunov ex- 
ponents for the system) is another chaotic property of 
many-particle swtems which has been studied exten- 
sively [H [S IB, H El, 0, 



21, 25, m, 27 



Each step in the exponent spectrum is associated with 
delocalised wavelike structures in the corresponding Lya- 
punov vector, now referred to as a Lyapunov mode. 

The significance of the step structure is that it appears 
in the Lyapunov exponents which are closest to zero, con- 
necting them with the slowest macroscopic modes of the 
system. Some analytical approaches, such as random 
matrix theory [soL 
riodic orbit theory 



311. kinetic theory [18|, |33| and pe- 
32l | have been used in an attempt to 
understand this phenomenon. A clue to understanding 
the stepwise structure of the Lyapunov spectrum and the 
Lyapunov modes is in the behavior of the Lyapunov vec- 
tors associated with the zero-Lyapunov exponents 21 1. 
In each case the Lyapunov modes associated with a zero 
exponent are Noether transformations [23| generated by 
the conserved quantities or time translation along the 
phase space trajectory. The Lyapunov modes associated 
with the steps in the spectrum are A:-vector analogues of 
the zero modes and have the same basis in the funda- 



mental symmetries of the system 1^ 2^ 21 1 . 

We consider a quasi-one-dimensional system (QOD) 
[lil where the two-dimensional rectangular system is 
so narrow that the particles remain ordered in the x- 
direction. We use Hard-wall boundary conditions in the 
X direction and Periodic boundary conditions in the y 
direction (the (H,P) boundary conditions). The step 
structure of the Lyapunov spectrum consists of one-point 
steps and two-point steps [l6|. The one-point steps cor- 
respond to transverse modes while the two-points steps 
correspond to longitudinal and momentum proportional 
modes. The significant advantage of the QOD system 
is that both the exponents and the modes can be ob- 
tained to high accuracy by standard numerical schemes 
0, H, 0, 0, H S 0] ■ For systems with smooth interac- 
tion potentials it has proved much more difficult to get 
clear numerical evidence for the steps in the Lyapunov 
spectrum and to find modes 13, 35 1. However, the same 
structure must exist as the dynamics is subject to the 
same invariances and conservation properties. 

The QOD system with (H,P) boundary conditions we 
studied contained TV hard disk particles. The complete 
description of the system at any time is contained in the 
4Af-dimensional phase vector {q,p) where q = qi,...qAr 
and p = pi,...pAr where and pi contain the x and 
y coordinates and momentum of particle i. The time 
evolution of the phase vector through a free flight, and 
then collision, is given by the matrix equation 
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where each scripted matrix is an N xN matrix containing 
2x2 sub matrices, 2 is the identity and O is the zero 
matrix. The matrix J\f that changes the momenta at 
collisions is given by 

// . . 0\ 
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Again, each element is itself a 2 x 2 sub-matrix. The 
term njj = {xij, yij) is a row vector containing the x and 
y components of the separation between particles i and 
j at collision, so that the dyadic product is a 2 x 2 

matrix. 

We write the Lyapunov vectors as {Sq, 5p)^ where 
5q and 5p are iV-dimensional vectors containing the 2- 
dimensional entries for each particle position separation 
5(\j or momentum separation (5pj. For this QOD sys- 
tem there are four zero-Lyapunov exponents and hence 
four associated Zero (Z) Lyapunov modes. The numeri- 
cally observed Z modes can be written as linear combi- 
nations of elements of the basis set [2lf] 5Ty = 0), 
5Vpy = (0,M, 5Vt = (p,0)/|b|| and 5T, = (0,p)/|b|| 
where the j*'' element of 5q and 5p in 5Ty and STpy 

is given by Sq^j ~ 5-pj ^ f ? V ''^hile p is an A^- 



1 



dimensional vector whose j*'* component is Pj (the mo- 
mentum of particle j) and ||p|| is the magnitude of the 
total momentum. 

The time evolution equations for tangent space dynam- 
ics consist of many repeats of a free-flight followed by a 
collision and the Gram-Schmidt procedure. The first two 
of these steps, the application of a free-flight matrix then 
a collision matrix, evolve the tangent vector in time from 
to T and can be written as 



5q{T) 
5p{t) 



where 
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Each component of the matrix Q is a 2 x 2 sub-matrix 
where the only non-trivial components are those associ- 
ated with the two particles that collide, i and j through 
Qij. Then Qij given is by 



Qt] = (n. 



I 



nij ■ pij 



P^jl 



n^j ■ Pi 



(5) 



where pfj = {pxij,Pyij) is a row vector containing the x 
and y components of the relative momenta at collision 
{Pij = Pj ^ Pi)- The principle property that we will 
exploit now is that Qij ■ = 0. 

To understand the tangent vector dynamics we con- 
sider the action of the matrix J\f on either Sq or Sp which 
gives 



AfSq = Sq 



((5qj - Sq,)X. 



(6) 



where X is the A^-dimensional column vector with all 
elements equal to zero except for Xi = 1 and Xj = —1. 



Note that for a QOD system j = 
result is general. 

Similarly, the action of the matrix Q qives 

QSq = Qij ■ {Sqj - Sqi)X 



1, but otherwise the 



(7) 



Any Lyapunov mode for which Sqi — Sqj = or Spi — 
Spj — exactly - such as the zero mode STy - is preserved 
by the dynamics, while the conjugate zero mode grows 
linearly with time, so for example STpy^r) = STpy + rSTy 

The application of the dynamics on the momentum 
dependent Z modes gives for STt that Afp = p and 
Qp = so that (p,0)/||p|| (p ,0)/\\p'\\ which is STt 
at the new time. Therefore the functional form of STt is 
preserved exactly. The conjugate mode is linear in time 
5Te{r) =ST, + TSTt. 

The numerically observed Lyapunov modes are of four 
types, Z modes, transverse (T) modes, longitudinal (L) 
modes and momentum proportional (P) modes. The L 
and P are usually observed together as a combined LP 
mode. The T modes are given by 



ST' = 



( Sq 
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where 7 and 7 are constants [23|. The components of 

both SqT and SpT are of the form Sqj = i ^ J where 

c-fij — COS k^iXj • 

In a transverse mode the term {Sc{j — Sqi) in equations 
([6]) and ([7|), becomes {Syj — Syi) = — e7„ sin(fc„Xi) where 
e = knXij = niTXij/L. Here Xij is the x-component of 
the distance between particles i and j at collision and L 
is the length of the system in the a;-direction. The small 
parameter e is linear in the mode number n. Thus the 
time evolution of a T modes is given by 



Sq{T) 
Sp{t) 



Sq + rSp 
Sp 



0{e). 



(9) 



and henceforth we neglect the e dependent terms. Clearly 
for sufficiently large e the dynamics without the order 
e term will become incorrect. We assume that below 
a threshold value of e the dynamics in equation ^ is 
correct. 

Further, for the LP mode, the P component is linear in 
the momentum, so that for example Sqi = PnPi 

COS kfiX'i 

then Sqi ~ Sq^ = /3„(p.y cos knX^ - epj sin knXi). As Qij ■ 
Pij — the first term is zero and the result is order e. 
The LP modes are of the form 
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where q;„ and are constants j23]. The A'^- 

dimensional vector pc has components PjC„j, Snt = 
sin LUnt and c„( = coscunt (with frequency w„). In the 
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case of (H,P) boundary conditions the LP exponents are 
doubly degenerate so there is a second LP mode i5LP^ 
with Snt and c„t interchanged which is orthogonal in time 
to SLP". This result means that the LP modes also have 
a time evolution governed by equation ^ with a differ- 
ent order e term and again we assume that this dynamics 
is correct below some threshold e. 

The numerically observed T modes are invariants of 
the dynamics and the pair of LP modes define a two- 
dimensional sub-space. The action of the dynamics and 
Gram-Schmidt procedure therefore simply results in the 
T mode remaining orthogonal to the centre space defined 
by the zero modes. This suggests that we can define an 
inside-out Gram-Schmidt procedure to obtain the same 
effect. To do this we make two assumptions: 1) that the 
functional forms for the T, L and P modes are known; 2) 
that the numerical values for the coefficients are known. 

Given these assumptions we can calculate the values 
of the Lyapunov exponents in the step region. If we con- 
sider a T or L mode then under a free flight and collision 
using equation ([9]) the mode changes, but then the Gram- 
Schmidt procedure returns it to its initial direction with 
some scaling factor (, as 



fSq 
[ 5p 
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The first right arrow is the action of the free flight and 
collision while the second right arrow is the result of 
the Gram-Schmidt procedure. Here we use a simplified 
inside-out Gram-Schmidt procedure in which we assume 
that the mode is already orthogonal to the centre space 
and then ensure orthogonality with respect to the conju- 
gate mode. 

Using the symplectic property of the system [2^ the 
mode conjugate to ST'^"'> is 



ST'- 



—Sp 
Sq 



(12) 



Applying the inside- out Gram-Schmidt procedure to 
equations (fTTj) gives 



c 



Sq 
Sp 



Sq + rSp 
Sp 



+ riSp-Sp)(^ fP'^ (13) 



or two equations to solve for the scale factor (. It is 
straigthforward to see that as the mode is normalised 
Sq ■ Sq -\- Sp ■ Sp — 1 so both components of equation (|13p 
give the same solution 



C = i 



_{Sp ■ Sp){Sq ■ Sq) 



(14) 



Sq ■ Sp 

The full time evolution is infinitely many repeats of this 
process: free-flight, collision and Gram-Schmidt, so the 
Lyapunov exponent is given by 



m 

A= hm T^lnJIC 



(15) 



TABLE I: A comparison of predicted (egns [T6] [TTl [T8)) and 
numerically observed Lyapunov exponents for a 200 particle 
QOD system with (H,P) boundary conditions at density p = 
0.8 and temperature T = \. We use units where the mass and 
disk radius K are 1, the total energy is A'^ and the system size 
is Ly = 1.157? and L^, = N / pLy where p is the density. There 
are small differences between the x and y projections of (3 so 
we include both results. 
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1 


0.0388 


0.0393 


0.0599 


0.0559 


0.0502 


0.0605 


2 


0.0749 


0.0784 


0.1169 


0.1044 


0.0965 


0.1229 


3 


0.1099 


0.1177 


0.1507 


0.1348 


0.1344 


0.1848 


4 


0.1431 


0.1571 


0.1753 


0.1472 


0.1591 


0.2484 


5 


0.1748 


0.1961 


0.1693 


0.1391 


0.1595 


0.3140 


6 


0.2007 0.2352 


0.1855 


0.1409 


0.1823 


0.3791 



assuming that 



For the mode (5T", Syj = jnCnj and Spyj 
,y ~ N/2, we have 

A™ = Y7„7n. 



7„c. 



(16) 



Similarly, we can consider the evolution of the nega- 
tive mode ST~^ and by ensuring orthogonality with 
respect to its conjugate mode (5T", the result will be 

A-„ = -^l-nl-n = — A„. 

Next treat the longitudinal part of the LP mode with- 
out its explicit time dependence. Clearly this is just the 
same as the transverse mode and the result is 



A. 



N 



(17) 



We treat the momentum dependent part of the LP 
mode without its explicit time dependence and ensure 
orthogonality with respect to the conjugate P compo- 
nent. A similar argument leads to the result 



A„ = iV/3„/3„r. 



(18) 

Ej Pj- Again 
'A„ for both L 



i=l 



where the temperature is given by 2A'^T 
the negative exponent is simply A_„ — 
and P components of the mode. 

In table U we compare the predicted and numerical re- 
sults for the Lyapunov exponents of the first six modes 
of each type. The first T mode and P mode are quite 
accurate but the first L mode is 20% less than the nu- 
merical result. Generally the results become worse for 
higher order modes. 

There are a number of possible sources of error. The 
dynamics is limited by the size of the neglected term e 
and will be worse for larger n. Probably the most impor- 
tant limitation is the simplified inside- out Gram-Schmidt 
procedure as it assumes that the functional forms for the 
modes are already orthogonal to the centre space of Z 
modes. At any finite N this is not correct and any more 
exact inside-out Gram-Schmidt procedure would need to 
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work systematically ensuring orthogonality with all pre- 
vious modes. Thus, for example, to Gram-Schmidt the 
3rd T mode it should be explicitly made orthogonal to 
the centre space, the 1st and 2nd T modes and any LP 
modes with lower value exponents. 

In conclusion, we have shown that the Lyapunov ex- 
ponents for all types of modes can be calculated using 
an inside- out Gram-Schmidt procedure and a complete 
knowledge of the functional form of the modes. The sim- 
plified Gram-Schmidt procedure is only accurate for the 
first T and P mode components but the systematic ap- 
proach suggested above may improve the accuracy but 
at the cost of the simplicity of the result. Thus we see 
that the same information that is encoded in the modes 
is also encoded in the values of the exponents. 

In all numerical calculations of Lyapunov modes there 
are a set of modes which are stable below some maximum 
mode number Uj-nax ■ Here we require that e in the dynam- 
ics of equation fO]) is less than some threshold Cmax which 
we can estimate from the numerics used to generate table 
m For a system of = 200 disks we find that Umax ^ 24, 
and for the QOD system Xij is positive and bounded by 

Wl - Lg/4 < < 1. Therefore ^max ~ O.SbXij which 

corresponds to about 8 particles per half wavelength. If 
the initial Lyapunov vector is mode-like for a particular 
n then the dynamics will preserve its mode-like charac- 
ter for n < Umax and it will be unstable for n > Umax- 
The question of the stability or otherwise of a particular 
mode is a different level of stability for this system. 
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